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Abstract 

A classical (non-quantum-mechanical) relativistic ideal gas in thermodynamic equi- 
librium in a uniformly accelerated frame of reference is studied using Gibbs's micro- 
canonical and grand canonical formulations of statistical mechanics. Using these meth- 
ods explicit expressions for the particle, energy and entropy density distributions are 
obtained, which are found to be in agreement with the well known results of the rel- 
ativistic formulation of Boltzmann's kinetic theory. Explicit expressions for the total 
entropy, total energy and rest mass of the gas are obtained. The position of the center 
of mass of the gas in equilibrium is found. The non-relativistic and ultrarelativistic ap- 
proximations are also considered. The phase space volume of the system is calculated 
explicitly in the ultrarelativistic approximation. 

PACS: 03.30+p, 05.20Gg, 05.20 Jj, 04.40-b 



Introduction: 

The relativistic ideal gas in a gravitational field has been extensively studied using the 
relativistic Boltzmann equation [H El El IH El El E] • This theory has found interesting appli- 
cations in astrophysics, cosmology and nuclear physics [6l [7] . 

General relativity is a local theory. The relativistic Boltzmann distribution function too 
is defined locally. It is therefore not surprising that Boltzmann's kinetic theory has been the 
method employed in the study of relativistic gases in gravitational fields. 

Gibbs's microcanonical formulation of statistical mechanics deals with systems of finite 
volumes and finite energies. To our knowledge, the methods of statistical mechanics, which 
are based on the microcanonical formulation, have not been employed so far to check the va- 
lidity of the results of the relativistic kinetic theory for a gas in the presence of a gravitational 
field. It seems of interest to check if both methods give the same results. 

In section 1 we study a classical relativistic ideal gas in thermodynamic equilibrium in 
a uniformly accelerated reference frame using Gibbs's microcanonical and grand canonical 
formulations of statistical mechanics. 

In section 2 we compare the results of section 1 with the well known results of the 
relativistic kinetic theory. We find complete agreement. 

In section 3 we show that the expressions for the particle, energy and entropy density 
distributions can be integrated and we obtain explicit expressions for the total entropy, 
total energy and rest mass of the gas. We also find the center of mass of the system in 
equilibrium and derive the condition for the relativistic gas to be in mechanical equilibrium 
in the uniformly accelerated frame. 

The non-relativist ic approximation is considered in the first part of section 4. The ap- 
proximate expressions for the total energy and total entropy are found to be in agreement 
with known results for a classical non-relativistic ideal gas in a uniform gravitational field. 

The ultrarelativistic approximation is studied in the second part of section 4. Explicit 
expressions are found in this case using the microcanonical formalism. In particular, the 
phase space volume of the system can be found explicitly in this case, which allows a direct 
comparison between the results of the microcanonical and the grand canonical formulations. 

Section 1: The microcanonical and grand canonical formulations 

Let us consider a gas consisting of a very large number N of identical (indistinguishable) 
structureless particles. Assume the particles of the gas do not interact with each other, except 
for elastic collisions among themselves and with the walls of the container (an ideal gas). 
Let us assume the gas to be adiabatically isolated: any changes in the total internal energy 
can only be the result of a change in the external parameters. We assume the container is 
in uniformly accelerated motion along the ^-direction with respect to an inertial reference 
frame K', such that the components of the four-vector acceleration x'^ (/i = 0, 1, 2, 3) of the 
container obey the condition: 

X' g = Tj^yX Q = ^ Oi (1) 

where is a constant, rj^y = diag{+l, —1, —1, —1) and c is the speed of light. 
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For simplicity, we assume the container to be a cyhnder of base area A and height L. 
The axis of the cyhnder is ahgned along the z-direction. 

In the local reference frame K of the accelerated observer, with respect to whom the box 
is at rest, the metric tensor can be written as [HJ [9] 0: 

(2) 

2 

We assume the cylinder containing the gas is in a region within a distance ^ of the 
observer. In other words, we assume the container is small enough and close enough to the 
origin of the reference frame K that the metric ^ is valid everywhere inside the container. 

The Lagrangian for a system of non-interacting identical relativistic particles in the 
non-inertial uniformly accelerated reference frame K can be written as: 

N 

L = -mc^i:((l + i5^) , (3) 

a=l 

where {i = 1,2,3) are the components of the velocity of the ath particle in K, = 
■^L + '^ay + ''^azy the mass of a particle. 

From the Lagrangian (jS]) we can find the components of the momentum of the ath particle: 

Pai = T, (4) 

1 + ^^' 

The Hamiltonian of the system can therefore be written as: 

where pI = pL + ply + pI,- 

Let us assume the gas has reached thermodynamic equilibrium. It is assumed that any 
infinitesimally small portion of the fluid contains a very large number of particles. 

In the non-inertial reference frame K, where all the portions of the fluid are at rest, the 
total energy of the fluid can be written as: 



E = aJ e{z)dz, (6) 



where e{z) can be interpreted as the energy density of the fluid measured with respect to K. 

In the instantaneous (freely falling) proper inertial frame K' , in which the fluid is also at 
rest, the total energy of the fluid can be written as follows: 



^In the non-relativistic limit this metric takes the form: g^^ = diag (l + —1, —1, — l), which is the 
expression for the metric describing a uniform gravitational field in this approximation [S]. 
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L L 

E' = aJ e'{z')dz' = aJ e'{z)dz, (7) 



where e' is the proper energy density (energy per unit proper volume) measured with respect 
to K'. Notice that both reference frames K' and K are at rest with respect to the fluid and 
with respect to each other, and the metric tensor in K is of the form ([2]), therefore dx' = dx, 
dy' = dy and dz' = dz. The dimensions of the container are the same in both frames of 
reference. 

In K' the particles of the gas move in rectilinear uniform motion in between collisions 
(as free particles), while in K they are subject to the effect of non-inertial forces. The two 
quantities e and e' are related as follows: 

l.f)/. (8) 

The energy density e of the gas (measured with respect to K) takes into account the 
effect of the non-inertial forces on the particles, while the proper energy density e' does not. 

In the non-inertial reference frame K, where all portions of the fluid are at rest, the 
components of the energy-momentum tensor of the fluid can be expressed as follows: 



\ 



T^"" = diag -2,P,p',P 



(9) 



where p' is the pressure in the fluid (a local function of position in the fluid) measured with 
respect to the instantaneous inertial rest frame K'. 

Equating the covariant derivative of the energy-momentum tensor to zero (law of energy- 
momentum conservation) we obtain the condition [101 [H] for the fluid to be in a state of 
mechanical equilibrium: 

dp' _ g {e'+p') ^^^^ 



dz c 



(l + f 



In the micro canonical formalism, the number Q{E, N, g, A, L) = D{E, N, g, A, L)5E of 
accessible microstates of the system (for a fixed number of particles and fixed external 
parameters g, A and L, and energy in the range {E, E + SE)) can be found from the phase 
space volume integral: 



$(i?, N,g,A,L) = j ■■■ j dxidyidzidpi^dpiydpiz---dxNdyNdzNdpNxdpNydpNz 

1 



a=l 

as follows: 
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The phase space volume integral f|TT]) is the volume in phase space enclosed by the surface: 

a=l _ 

In ( [T2|) . h is Planck's constant and the A^! accounts for the fact that the particles of the 
gas are indistinguishable. 

From (ITT]) we immediately find: 

^{E,N,g,A,L) = {AnAf J ... J pl..p%dzi...dzNdpi...dpN 

l(,^E.)(r,jLf<^. (13) 



a=l 



Assuming that the number of particles is very large (A^! ^ y2^(f) (Stirling's 
approximation)), the entropy of the isolated system in thermodynamic equilibrium can be 
found from ffT2l) and ffT3l) as: 



S = k\nn{E,N,g,A,L), (14) 

where k is Boltzmann's constant. 

The quantities T and ;U can be found from the partial derivatives of the entropy as follows: 



f ~ [dE 



(15) 

N,g,A,L 

[dN 



I ■ (16) 

E,g,A,L 

Let us consider a layer of the fluid (of area ^4) at a height z. This is an open subsystem. 
For this portion of the gas we can write the grand canonical partition function: 



iV „ / 2 \i 

^ H-(Adz)^ f f -Ef^(i+f)(i+J^) . . 



7V=0 



(oo \ 
----(x)7-"'^'"-''*H 

exp I An Adz i^!!^^ J e-''^^^('^'")K'dx\ , (17) 
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where: 



m(? ( gz 



kT \ 6 

We notice that: 



«W = T7f(1 + %)- (18) 



1 e-<^+"')'a;2rfx = j u^^^^e-'^'^ du = (19) 

1 ^ 

where K2 is the modified Bessel function [121 [13] . 

From firT]) . flTK]) and f[T^ we finally obtain the grand canonical partition function for the 
relativistic ideal gas: 



fmcY K,{^ (l + g)) 



Z = exp e^AivAdz — ^ . (20) 



\ h J 



The grand canonical potential per unit volume in K, can be found directly from ( |20|) : 

, /mcX3^^{#{l+jf)) 

w = -kekTin — \ / r^- (21) 

ff^(l + f) 

From ( |2Tll the entropy density s (entropy per unit volume) and the particle density n 
(number of particles per unit volume) can be found using the equations: 



(22) 



dT ' 



n = -|'^) ■ (23) 

9,T 



Using the thermodynamic identity: 



u = e — Ts — fm, (24) 

and expressing fi in terms of the particle density n(0) at the bottom of the container, we 
obtain from (1^ . 023), ^^3^ and the particle density, entropy density, and energy density 
equilibrium distributions: 

n(.) = .(0) f^^('+/)| , (25) 

.(.)^.(.). in|i^(^r^(^(i±# 

'n(z)V/iy #(i + f) y V ir2(f#(i + f)) 

(26) 
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.(.)^ni.)kTU-'^(l^?i)'im^^], (27) 




where < 2: < L (we assume the origin of the reference frame K to be at the bottom of the 
container). 

The pressure in K can be found directly from the expression of the grand canonical 
potential per unit volume as follows: 

p = -u = n{z)kT (28) 
From fl27|) and ([8]) we find the proper energy density as a function of z: 

Substituting fl29l) into ffTOl) . we can solve the differential equation for the pressure p' in 
K' to obtain: 

(1 + S) 

The thermodynamic relations between e', ra', s', p' and the local temperature T' measured 
in the instantaneous inertial frame K' for a layer of gas in equilibrium at height z can be 
found from the grand canonical partition function (in the inertial frame): 



exp|e5*ft4x^<(4!^)'^iMk|. (31) 



As before, we can find s' and n' from the equations: 



'■-iw^ . (32) 



where u' can be found from fIST]) : 



= -fce#^4. ( 7: ) ^A^. (34) 

feT' 
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Finally, using the thermodynamic identity: 



u' = ^ - T's' - ii'n', 



and the relation: 



u' = -p 



(35) 



(36) 



for a relativistic ideal gas (in the inertial frame K'), in the absence of external forces, one 
finds the Juttner relations [H HSl [ISl El HSl US [20] : 



t' = p \ l 



mc 



p' = n'kT 



nk 



In 



n' \ n / I kV K2 (0) 

Comparing ([37D, ([38D and dSQ]) with ([29D, ([30]) and ([26]) we obtain: 

T 



T'{z) 



n = n, 



(37) 
(38) 

(39) 



(40) 
(41) 



s' = s. (42) 

Equations (14T!) and fjl2]) can also be obtained by noticing that the number of particles 
and the entropy are invariant quantities (have the same values both in K and K') and the 
volumes in K and K' are the same (as noted above). 

Equation fj^Ol) is the well known Tolman relation [2T|, [9] {T'^/g^ = T = const) between 
the local temperature T'{z) measured with respect to the inertial frame K' and the constant 
T defined in the non-inertial frame K. We have arrived at ( HO]) using the methods of 
statistical mechanics. 

Section 2: The relativistic kinetic theory 

One can also obtain Eqs. ([25] - [27]) from the kinetic theory, using the relativistic Boltz- 
mann equation [D [2] [S] [2 [S] [Q [7] . 

In the non-inertial uniformly accelerated reference frame K, the relativistic Boltzmann 
equation can be written as follows: 
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dt ^^^^^V\ Of J \ m^cy dp. 



(1 + 



— J u {f;p,pi;p,p^^ (/ {t,f,pj f (t,f,p^^ - f {t,f,p) f {t,f,pi)^ d'^pid^pd^pi. 

(43) 



The right hand side of the relativistic Boltzmann equation (H3l) is the colhsion integral 
It gives the rate of change, due to colhsions, in the mean number of particles in the phase 
space volume dxdydzdpxdpydp,. The left hand side of P3|) is ^ = ^ + {/, H}, where {, } 
is a Poisson bracket. The particular expression in the left hand side of fj43l) is obtained by 
inserting the Hamiltonian ([5]) in the Poisson bracket. 

For elastic collisions: 

P + Pi = P + Pi, (44) 

u,j^)%(,^jn'j,^4\\h^4\\ (45) 

The mean number of particles per unit volume is given by the equation: 

n(t,r) = -^ / /(^'^'^^^V (46) 

In equilibrium: 



dt 



(47) 



/eg (r, p) feq [t, Pi) = feq (r, p) feq {v, Pl) . (48) 

In Eq (HHj) it is assumed that the initial and final momenta obey the conditions (jHl HS]). 
In our case, due to the symmetry of the problem considered, we also have the following 
conditions on feq-. 



dfeq dfeq 



0. (49) 



dx dy 

From ( l43|) . (!47|-|49]) we obtain that the distribution function feq satisfies the equation: 



m 



(1 + 



dz \ m?c^ J dpz 



The solution of ( 150|) . which also satisfies the condition (jUj) (assuming (H ^ H5i) ). can be 
written as: 
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1 



feq{z,p} = ae " 'V - c y ^ (^5;L) 

where a and /3 are constants. 

Substituting ( |5T]) into ( H6l) and performing the integration we obtain: 



n{z) = ^ ^ . " . 52 

^ fl + !f 



In terms of the particle density at the bottom of the container we obtain: 

„(0) A',(/jmc'(l + g)) 



The energy density e{z) in can be derived from the relation: 



1 



e{z) = ljrnc^{l + f ) + /e. i^^P) d'p (54) 

Substituting fl5T|) into fl5^ . taking into account fl52|) . and performing the integration we 
obtain: 

,M "Wfi ,V.,^ri|^-n ^^^^(^""Hl + g)) ^ (55) 



Finally, we can obtain the entropy density from the relation |22] : 



(56) 



Substituting flSTl) into fl56l) . taking into account fl52|) . and performing the integration we 
obtain: 



s{z) = n{z)k 



V / /3mc2 (l + f ) 



/3mc^ 1 + ^ 



Comparing ([531 [541 EZD with ([251- [27D we find the constant /3: 



K-(/?mc^ (l + g)) - 

c'J K2{^mc^ (l + f)) 
(57) 



(58) 



Section 3: Global quantities 

The density of particles at the bottom of the container n(0) can be related to the total 
number of particles in the system by direct integration of ( !25l) : 



9 



(59) 

where we used the relation: ^ = -£ (^) 

The total energy of the system E in K can be found by substituting ( 125|) and (l59i) into 
(p7|) and performing the integration ([6]): 



= NkT + Nmc 




Similarly, substituting (|25|) and f l59|) into fl26|) . and integrating, we obtain the total en- 
tropy of the gas: 



S . Mk^M {^\' f^V^ U - "'i^^i'^.^)) 



N \ h J \ mc^ 



W {K. (#) - (l + f 



/ 2N ■^lfT#(l+4)' 



(61) 



The total energy of the system E' in K' can be found by substituting fl2S]) and flS^ into 
and performing the integration ([7j): 



i^.(#(i+f)) ' 

(1+4) 

= A^mc^^ 7 ^ r^- (62) 

^T#(i+f))' 



ii^' given by Eq (l62|) is the proper energy of the gas. 

The position of the center of mass of the gas with respect to the inertial frame K' can 
be determined by the equation [23]: 

L 

A f ze'{z)dz 

The integration in ( 163|) can be easily carried out and using (l62l) and ( 160|) one finds that: 

Mgz,,^, = E-E'. (64) 
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where the rest mass M of the fluid is defined as: 



M='^ (65) 
From (ES]), (1621), (EHl), (EHl) and (159^ notice that: 

Mg = A{p{Q)-p{L)). (66) 

The above equation can be interpreted as the condition (in K) for the object to be in 
static mechanical equilibrium: the net force exerted by the walls of the container on the gas 
as a whole equals its weight. Notice that in the non-inertial reference frame K the net force 
exerted on the gas by the side walls of the cylinder is equal to zero. 

Section 4: The non-relativistic and ultrarelativistic approximations 

The non-relativistic approximation. 

In the non-relativistic limit — > oo) we can use the expressions for the asymptotics 
of the modified Bessel functions 



for = 1,2 in fl60|) and fl6T]) to obtain the non-relativistic expressions for the total energy 
and total entropy of the gas: 

E = Nm^ + \NkT- /!^'^ (68) 
^ (e fcT — Ij 



mgL 1 \ .ri ^ ( A{2nmkT)l ^kT 



mgL ~ 



S = Nk\l- l^p^ j + «>n y-^-^e.- (1 - e-- ) I . (69) 

Equations fl68l) and fl69|) are in perfect agreement with the results for the non-relativistic 
classical ideal gas in a uniform gravitational field obtained in |24j. 

In the non-relativistic limit, where (1 + ^^^^j ^ ~ 1 + ^^2^2 1 the integration in f lT5]) can 
also be carried out and (assuming E > Nmc^ + NmgL) we obtain results in agreement with 
[25]. 

The ultrarelativistic approximation. 

In the ultrarelativistic approximation, where (1 + ■:;^2^] ~ the integration in (1T3|1 
can be carried out explicitly and (assuming E > Nmc^ (^1 + ^^) we obtain: 



^E,N,g,A,L) = -—\^E' 
(3iV)! go 



^ ^ N 

AirA 



( 



2 



(70) 
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From fITOl) . f|T2l) and f|T^ . using Stirling's approximation, we obtain a formula for the 
entropy of an ultrarelativistic ideal gas in thermodynamic equilibrium in a uniformly accel- 
erated reference frame: 



In the ultrarelativistic limit — 0) we 
for the modified Bessel functions as 2: — )• 0: 




) 



u 



(72) 



2 \z 



for z/ = 1, 2, in (160|) and ( 161]) to obtain: 



(73) 



47re^ A 




^ = A^A; In 



(A;T)^ 1 



(74) 



V 



Combining Eqs. ( !74|) and (!73|) we obtain (!7T]) . which shows that the microcanonical and 
the grand canonical formalisms give identical results in the ultrarelativistic approximation. 

In this paper we have considered the case of an ideal gas in a uniformly accelerated frame. 
The relativistic kinetic theory is applicable in the more general case of a gas in a stationary 
gravitational field. It would be of interest to generalize the results obtained here to find the 
microcanonical and grand canonical formulations for the gas in a stationary gravitational 
field. 
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